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■Abstract 

Event  rates  encountered  by  satellite-borne  sensors  typically 
exceed  the  available  telemetry  bandwidth,  necessitating  on¬ 
board  data  analysis.  One  method  of  reducing  output  telemetry 
requirements,  employed  by  data  processing  units  (DPUs)  on 
several  recent  and  forthcoming  satellites  (VIKING,  CRRES, 
CAMMICE)  is  to  partition  the  underlying  event  parameter 
spaces  into  disjoint  classes,  and  to  downlink  statistics  based  on 
these  classes  rather  than  individual  events.  The  DPUs  utilize  a 
hardware  implemented  binary  search  algorithm  to  classify  each 
valid  event  into  one  of  several  mass  and  mass  per  charge-state 
(M/Q)  groups,  using  look-up  tables  written  under  software 
control.  Mass  and  M/Q  group  boundaries  are  defined  via 
polynomial  curves  in  the  underlying  parameter  spaces.  The 
boundary  curves  were  obtained  using  a  combination  of  global 
and  local  least-squares  approximation  algorithms  applied  to 
pre-flight  sensor  calibration  data. 

I.  Introduction 

Event  rates  encountered  by  satellite-borne  sensors  typically 
exceed  the  available  telemetry  bandwidth,  necessitating  on¬ 
board  data  analysis  and  compression.  One  method  of  reducing 
telemetry  requirements  is  to  compute  summary  statistics  on¬ 
board  and  to  downlink  these  statistics  rather  than  individual 
events. 

For  example,  the  microprocessor-based  data  processing 
units  (DPUs)  that  support  the  Magnetospheric  Ion  Compo¬ 
sition  Spectrometer  (MICS)  aboard  both  the  Swedish  satellite 
VIKING  and  the  Combined  Release  and  Radiation  Effects 
Satellite  (CRRES)  utilize  a  hardware  implemented  binary 
search  strategy  to  classify  each  valid  event  into  one  of  several 
“mass  groups”  based  on  energy  (E)  and  time-of-flight  (T) 
values  generated  by  the  sensor.  Each  event  is  similarly 
assigned  to  a  mass  per  charge-state  group  based  on  the  time-of- 
flight  and  energy  per  charge-state  (E/Q)  of  the  incident  particle. 
The  assigned  mass  (M)  and  mass  per  charge-state  (M/Q) 
groups  are  then  used  to  index  a  matrix  of  hardware  accum¬ 
ulators  (called  “scalers”)  in  such  a  way  that  each  event  is 
counted  by  exacuy  one  M  vs.  M/Q  “matrix”  scaler.  By 
periodically  downlinking  the  complete  set  of  M  vs.  M/Q 
matrix  scalers,  a  coarse  determination  of  event  populations  can 
be  obtained  through  the  limited  bandwidth  telemetry  channel. 

Finer  temporal  and  parametric  resolution  is  achieved 
through  the  use  of  individual  “rate  scalers"  which  may  be 
programmed  to  accept  only  those  events  having  particular 
combinations  of  M  and  M/Q  values.  The  additional  inclusion 
in  the  telemetry  output  of  a  small  number  of  “direct  events” 
(raw  energy  and  time-of-flight  values)  allows  resolution  at  the 
scale  of  a  single  event  The  combined  use  of  M  vs.  M/Q 
matrix  scalers,  rate  scalers,  and  sample  direct  events  thereby 
provides  coverage  of  event  parameter  space  at  a  range  of  reso¬ 
lutions,  from  coarse  to  fine. 


The  following  sections  describe  in  detail  the  data  analysis 
and  reduction  schemes  employed  by  the  data  processing  unit 
(DPUA)  supporting  the  MICS  sensor  aboard  the  CRRES 
satellite.  Similar  techniques  are  utilized  by  the  DPUs 
associated  with  the  VIKING  MICS  instrument  [1]  and  the  Low 
Energy  Magnetospheric  Ion  Composition  Sensor  aboard 
CRRES.  Data  from  the  MICS  sensor  to  be  flown  on  the 
forthcoming  CAMMICE  satellite  will  also  be  analyzed  via  M 
vs.  M/Q  reduction. 

Classification  of  E-T  events  into  mass  and  M/Q  groups  is 
performed  automatically  by  a  hardware-implemented  binary 
search  algorithm  utilizing  look-up  tables  written  by  the  DPU 
software.  The  tables  contain  discrete  approximations  to  poly¬ 
nomial  curves  which  separate  adjacent  groups.  The  curves  were 
derived  from  extensive  pre-flight  sensor  calibration  data  using  a 
novel  curve-fitting  and  extrapolation  algorithm  combining 
local  and  global  least-squares  approximations.  Both  the  on¬ 
board  table  writing  software  and  the  techniques  employed  to 
derive  the  boundary  curves  are  described  in  detail  below. 

II.  DPU  Hardware 

The  data  processing  unit  (DPUA)  that  interfaces  with  the 
CRRES  MICS  instrument  performs  the  following  functions: 
(1)  Accepts,  decodes,  and  executes  ground  commands  uplinked 
to  the  DPU  by  way  of  the  spacecraft;  (2)  Distributes  power  to 
the  sensor  and  controls  the  high  voltage  levels  of  the  MICS 
electrostatic  analyzer  and  post-accelerator;  (3)  Controls  the 
experiment  by  configuring  the  sensor  in  particular  data 
collection  modes  through  the  use  of  “MICS  instructions”;  and 
(4)  Reads,  processes,  and  outputs  to  telemetry  data  from  the 
MICS  instrument.  The  following  sections  are  concerned 
primarily  with  the  last  of  these  functions;  additional  infor¬ 
mation  on  the  remaining  functions  may  be(  found  in  [2], 

A.  Overview'  of  PPUA 

As  shown  in  Fig.  1,  DPUA  consists  of  a  computer,  peripheral 
sections  such  as  various  I/O  buffers  and  the  MICS  Event 
Processor,  a  power  supply,  and  an  analog  signal  handler.  The 
computer  consists  of  a  microprocessor  (SANDIA  1802), 
control  logic,  and  memory  (12K  bytes  of  ROM  and  16K  bytes 
of  RAM).  The  radiation-hardened  1802  microprocessor  operates 
with  a  1.6  MHz  clock.  The  DPUA  software  is  initially  stored 
in  ROM  and  is  copied  to  RAM  shortly  after  the  DPU  is  turned 
on.  Power  to  the  ROMs  is  then  shut  off,  and  thereafter  the 
software  executes  exclusively  in  RAM. 

One  important  group  of  peripheral  sections  functions  as  the 
MICS  Event  Processor  (see  the  area  inside  the  dashed  lines  in 
Fig.  1)  which  analyzes  MICS  sensor  data  semi-autonomously 
after  initial  preparation  (table  writing  by  the  computer)  has 
been  performed.  The  MICS  sensor  generates  energy  (E)  and 
time-of-flight  (T)  values  for  each  event.  In  the  Event  Processor 
the  E  and  T  values,  along  with  an  energy  per  charge-state 
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DPUA  BLOCK  DIAGRAM 


Fig.  1.  DPUA  Functional  Block  Diagram. 


(E/Q)  value,  are  used  to  classify  the  event  into  one  of  several 
mass  (M)  and  mass  per  charge-state  (M/Q)  groups.  This 
classification  is  performed  by  the  hardware-based  Successive 
Approximation  Routine,  which  utilizes  look-up  tables  written 
by  the  DPU  under  software  control. 

The  assigned  M  and  M/Q  values  for  an  event  are  used  to 
address  one  scaler  in  a  5 12-element  “matrix”  of  scalers  (M  vs 
M/Q  in  Fig.  1).  DPUA  also  contains  eight  "rate  scalers"  in 
which  events  with  special  combinations  of  M  and  M/Q  are 
accumulated.  The  Rate  Scaler  Select  Table  determines  which 
combinations  of  M  and  M/Q  are  accumulated  by  the  individual 
rate  scalers.  Sample  "direct  events"  (E  and  associated  T  values) 
are  collected  from  the  sensor  and  stored  in  the  event  buffer. 
Direct  events  are  selected  according  to  a  priority  scheme  deter¬ 
mined  using  DPUA’s  Event  Priority  Table.  The  contents  of 
the  various  scalers  and  the  direct  event  buffer  are  read-out 
periodically  by  the  computer,  processed  (e.g.,  compressed),  and 
placed  in  the  telemetry  stream. 


For  each  valid  event,  the  MICS  sensor  outputs  an  8-bit  energy 
(E)  value  and  a  10-bit  time-of-flight  (T)  value.  These  values, 
together  with  a  S-bit  quantization  (E/Q)  of  the  high  voltage 
level  of  the  MICS  electrostatic  analyzer,  are  used  by  DPUA’s 
MICS  Event  Processor  to  classify  the  event  into  one  of  up  to 
16  mass  groups  (M)  and  one  of  up  to  32  mass  per  charge-state 
groups  (M/Q).  Mass  and  M/Q  group  classification  is  per¬ 
formed  directly  by  the  DPU  hardware  using  look-up  tables 
written  under  software  control.  Each  mass  group  is  defined  as 
the  region  of  E  vs.  T  parameter  space  lying  between  two 
predefined  boundary  polynomial  curves  of  the  form  T  =  P(E), 
as  discussed  further  in  Section  IV.  Soon  after  the  DPU  is 
powered-up  in  space,  the  DPU  software  calculates  discrete 
approximations  to  these  curves,  based  on  stored  coefficients, 


and  writes  the  calculated  values  to  an  area  of  DPU  memory 
known  as  the  Mass  Croup  Table. 

The  Mass  Group  Table,  or  MGT,  is  a  two-dimensional 
array  indexed  by  mass  group  M  (4  bits),  and  energy  E  (8  bits). 
For  each  mass  group  M,  the  256  table  values  MGT(E=0,M) 
...MGT(E=255,M)  give  the  discrete  approximation  to  the 
upper  boundary  polynomial  for  the  mass  group  -  that  is, 
MGT(E,M)  is  the  nominal  time-of-flight  value  (10  bits)  for  a 
particle  of  energy  E  belonging  to  mass  group  M.  The  DPU’s 
hardware-based  Successive  Approximation  Routine  (SAR)  uses 
the  precomputed  Mass  Group  Table  values  to  classify  each 
incoming  event  (specified  by  E  and  T)  via  a  binary  search  for 
the  mass  group  M  satisfying  the  inequality  MGT(E,M)  £  T  < 
MGT(EM+1)-  In  particular,  the  SAR  first  compares  T  to  the 
table  value  MGT(EJM=8),  i.e.,  to  the  table  value  corre¬ 
sponding  to  E  for  mass  group  number  8.  Depending  on 
whether  T  <  MGT(E,M=8)  or  T  2  MGT(E,M=8).  the  SAR 
then  compares  T  to  MGT(E,M=4)  or  MGT(E>f=12),  respec¬ 
tively.  The  first  comparison  determines  in  which  half  of 
“mass  group  space”  the  event  belongs;  subsequent  compar¬ 
isons  (to  MGT(E,M=4)  or  MGT(E,M=12),  for  example) 
determine  the  appropriate  quadrant,  octant,  and  sextant.  Thus, 
in  at  most  four  comparisons,  the  unique  mass  group  M  is 
found  which  satisfies  MGT(E,M)  ST<  MGT(EJvf+l);  the 
SAR  then  assigns  the  event  being  classified  to  mass  group  M, 
as  illustrated  in  Fig.  Z 

While  this  classification  algorithm  imposes  slighdy  greater 
hardware  processing  demands  than  would  direct  table  look-up 
(storing  the  mass  group  directly  as  a  function  of  E  and  T),  the 
reduction  in  memory  requirements  is  quite  dramatic:  (16 
groups)  •  (256  E-values/group)  •  (10  bits/T-vaiue)  =  40,960 
bits  for  the  SAR-based  scheme,  versus  (256  E-values)  *  (1024 
T- values)  •  (4  bits/group)  =  1,048,576  bits  for  direct  table 
look-up.  Thus  the  slight  increase  in  hardware  complexity  is 
more  than  offset  by  the  96%  decrease  in  the  amount  of 
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memory  required  (it  should  be  noted  that  the  throughput  of  the 
SAR  exceeds  the  event  rate  of  the  M3CS  instrument,  hence  no 
bottleneck  is  introduced  by  this  increase  in  processing). 

The  M/Q  Group  Table  is  organized  in  a  similar  fashion. 
The  values  in  this  table  represent  discrete  approximations  to 
polynomial  curves  of  the  form  T  =  P(E/Q),  where  E/Q  is  a 
quantization  of  the  energy  per  charge-state  of  the  event.  Up  to 
32  M/Q  groups  can  be  defined,  so  that  the  table  requires  a  total 
of  32  •  32  •  10  =  10,240  bits,  as  opposed  to  32  •  1024  •  5  = 
163,840  bits  for  a  direct  table  look-up  approach  -  the  reduction 
in  storage  requirements  is  thus  approximately  94%. 

III.  On-Board  Data  Processing  Schemes 

The  principal  sensor-related  data  entries  processed  by  the  DPU 
are  the  R0-R7  rate  scalers  and  the  M  vs.  M/Q  matrix  scalers. 
The  contents  of  scalers  (RAM  accumulators)  are  defined  via 
software-configurable  look-up  tables  residing  in  DPU  memory. 
The  DPU  software  periodically  reads  and  resets  the  scalers, 
compresses  the  values  read  from  up  to  19  significant  bits  to  8 
bits,  and  then  places  the  compressed  data  in  the  telemetry 
stream.  The  matrix  and  rate  scalers  are  described  in  turn 
below;  a  detailed  description  of  the  scaler  compression  scheme 
utilized  by  DPUA  may  be  found  in  [2]. 

A.  M  vs.  M/Q  Matrix  Scalers 

As  discussed  above,  each  valid  sensor  event  is  assigned  to  one 
of  up  to  16  mass  (M)  groups  and  one  of  up  to  32  mass  per 
charge-state  (M/Q)  groups  by  DPUA’s  MICS  Event  Processor, 
based  on  the  contents  of  the  Mass  Group  and  M/Q  Group 
Tables.  The  DPUA  hardware  maintains  a  16  x  32  =  512- 
element  matrix  of  scalers  (RAM  accumulators),  indexed  by  M 
and  M/Q,  respectively.  A  given  M  vs.  M/Q  matrix  scaler 
counts  only  those  events  assigned  to  the  corresponding  mass 
and  mass  per  charge-state  groups.  The  scalers  are  19  bits  each, 
allowing  up  to  2l$-l  events  to  be  accumulated  without  over¬ 
flowing  (in  practice,  these  scalers  never  overflow).  Prior  to 
telemetry  output,  all  scalers  are  compressed  to  8  bits. 

Notice  that  the  mapping  from  energy  and  time-of-flight 
values  (sensor  outputs)  to  M  and  M/Q  groups  (Event  Proces¬ 
sor  outputs)  effectively  reduces  the  dimensionality  of  event 
parameter  space  from  8  •  1024  to  16  •  32,  or  approximately 
94%.  By  periodically  downlinking  the  entire  complement  of 
M  vs.  M/Q  matrix  scalers,  a  coarse  determination  of  event 
populations  can  be  obtained  through  the  limited  bandwidth 
telemetry  channel. 

B.  Rate  Sealers 

The  outputs  of  the  Event  Processor  (M  and  M/Q)  can  also  be 
stored  in  eight  rate  scalers  that  are  output  to  telemetry  much 
more  frequently  than  the  M  vs  M/Q  matrix  scalers.  Each  rate 
scaler  (R0  -  R7)  may  be  “programmed”  to  accept  only  those 
events  with  particular  combinations  of  M  and  M/Q.  This 
scheme  allows  multiple  values  of  M  and  M/Q  to  be  grouped 
together  and  counted  by  a  single  scaler. 

The  assignment  of  (M.M/Q)  combinations  to  particular 
scalers  is  accomplished  using  the  software -configurable  Rate 
Scaler  Select  Table.  There  are  eight  distinct  “pages"  in  the 
Rate  Scaler  Select  Table,  each  of  which  defines  a  different 
assignment  rule  for  the  eight  rate  scalers.  Only  one  page  at  a 
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Tn  =  TABLE  VALUE  *  1024  -  Tn 
T  2  Tn  <S>T  +  !„  >  1023 


Fig.  2.  SAR  Mass  Group  Classification:  Incoming  events 
specified  by  energy  (E)  and  time-of-flight  (T)  are  classified  in 
one  of  sixteen  mass  groups  (M)  by  successive  comparisons  of 
T  with  the  nominal  time-of-flight  values  stored  in  the  DPU’s 
Mass  Group  Table.  Shown  is  a  schematic  representation  of 
table  values  for  a  particular  value  of  E. 


time  is  active,  i.e.,  defines  the  contents  of  the  rate  scalers  (the 
currently  active  page  may  be  changed  by  uplinking  the 
appropriate  ground  command  to  the  DPU).  The  rate  scalers 
have  a  capacity  of  2I6-1  (16  bits)  per  scaler  and  are  compressed 
to  8  bits  prior  to  telemetry. 

Events  arriving  at  the  Event  Processor  are  automatically 
routed  to  the  appropriate  rate  scaler  based  on  the  definitions  of 
R0  -  R7  contained  in  the  currently  active  page  of  the  Rate 
Scaler  Select  Table.  Each  of  the  eight  pages  can  be  viewed  as 
a  matrix  indexed  by  M  (4  bits)  and  M/Q  (5  bits).  Each  location 
within  this  32  x  16  matrix  contains  a  3-bit  value  in  the  range 
0  to  7  specifying  the  rate  scaler  (one  of  R0  -  R7)  in  which 
events  with  the  corresponding  M  and  M/Q  values  are  to  be 
accumulated. 

In  this  way,  regions  of  particular  interest  in  M  vs.  M/Q 
space  can  be  isolated  and  monitored  with  greater  parametric  and 
temporal  resolution  than  is  possible  using  the  matrix  scalers 
(the  accumulation  period  for  the  DPUA  rate  scalers  is  1 /256th 
that  of  the  matrix  scalers). 
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TV.  Mass  and  MIQ  Group  Determination 

As  discussed  above,  DPUA’s  MICS  Event  Processor  classifies 
each  sensor  event  into  the  mass  CM)  and  mass  per  charge-state 
(M/Q)  groups  determined  by  values  contained  in  the  Mass  and 
M/Q  Group  Tables,  respectively.  The  contents  of  the  Mass 
Group  Table  (for  example)  represent  discrete  approximations  to 
polynomial  curves  of  the  form  T  »  P(E),  where  T  and  E  are 
nominal  time-of-flight  and  energy  values.  These  curves  divide 
the  underlying  E  vs.  T  parameter  space  into  disjoint  regions, 
each  of  which  constitutes  a  single  mass  group. 

In  this  section  the  genesis  of  die  mass  group  boundary 
polynomials  is  described;  similar  comments  apply  to  the  M/Q 
group  boundary  polynomials.  An  overview  of  the  techniques 
used  to  obtain  the  polynomial  coefficients  is  presented  first, 
followed  by  a  brief  description  of  FIT,  which  is  an  integrated 
software  package  implementing  the  procedures  outlined  in  the 
overview.  Finally,  the  mathematical  foundations  of  the  curve¬ 
fitting  algorithms  are  reviewed. 

A.  Overview 

The  boundary  curves  separating  adjacent  mass  groups  are 
derived  from  polynomial  least-squares  fits  to  calibration  data 
obtained  from  the  MICS  instrument  prior  to  launch.  At  first 
glance  it  might  appear  that  one  need  only  calculate  a  single 
approximating  curve  for  each  ion  species  and  then  translate 
these  curves  slightly  so  that  the  curves  separate,  rather  than  fit, 
the  data  corresponding  to  distinct  ion  species.  However,  the 
curves  so  obtained  tend  to  intersect  for  low  energy  values  and 
are  therefore  unsuitable  for  use  with  an  S AR-based 
classification  scheme.  In  particular,  for  each  pair  of  E  and  T 
values,  the  SAR  requires  that  there  be  a  unique  mass  group  M 
satisfying  the  inequality,  MGT(E,M)  ST<  MGT(E,M+1), 
which  is  equivalent  to  the  requirement  that  the  mass  group 
boundary  curves  be  non-intersecting.  In  order  to  satisfy  this 
constraint,  a  model-based  approach  was  employed  to  determine 
the  polynomials. 

Assuming  mass,  energy,  and  time-of-flight  to  be  related  by 
a  multivariate  polynomial  of  the  form  M  =  Q(E,T)  (e.g. 
Q(E,T)  ®  Po(E)  +  T*Pj(E)  +  T2-P2(E),  where  Po,  Pi,  and  P2 
are  polynomials  in  E),  a  “global"  least-squares  approximation 
was  performed  over  die  entire  calibration  data  set,  resulting  in 
the  determination  of  the  polynomial  coefficients  relating  E  ami 
T  to  the  atomic  mass  M  of  each  ion  species. 

The  relationship  summarized  by  M  =  Q(E,T)  was  then 
used  to  interpolate  data  points  for  “fictitious”  values  of  M  by 
letting  E  range  from  0  to  255  and  solving  M  =  Q(E,T)  for  the 
corresponding  values  of  T.  For  example,  with  Q(E,T) 
above,  evaluation  of  Po,  Pi,  and  P2  at  a  particular  value  E  = 
E'  yields  three  values:  Eo  =  Po(E),  Ei  =  PifE1),  and  E2  - 
P2(E’).  Substituting  these  values  into  the  equation  for  M  then 
yields:  M  =  Eo  +  Ei*T+  E2*T2,  or  equivalently,  0  =  (Eo  -  M) 
+  Ei*T  +  E2*T2,  where  the  Et’s  are  constants.  After  setting 
M  to  some  “fictitious"  value  (e.g.,  M  =  150,  the  quadratic 
equation  0  =  (Eo-M)  +  Ei*T  +  E2*T2  can  be  solved  for  T. 

By  fixing  M  and  letting  E  vary  from  0  to  255,  a  data  set  of 
“simulated"  E-T  pairs  satisfying  M  =  Q(E,T)  was  obtained,  as 
described  above.  A  polynomial  least-squares  fit  to  this  data  set 
was  then  performed,  producing  a  curve  T  =  P(E)  ®  ao  +  ai*E  + 
a2*E2  +  a3*E3,  describing  the  nominal  relationship  between  E 


Ion  Energy  (chonnels) 

Fig.  3.  Mass  Group  Boundary  Curves  Calculated  by  FIT. 

and  T  for  the  “fictitious”  mass  value  M.  (In  practice  the 
logarithms  of  E  and  T  were  utilized.)  For  example,  with  M  = 
1.5,  a  curve  separating  the  data  for  M  =  1.0  (hydrogen)  from 
the  data  for  M  =  2.0  (deuterium)  was  obtained.  This  technique 
was  repealed  for  various  "fictitious”  mass  values,  producing  a 
set  of  non-intersecting  mass  group  boundary  curves  suitable 
for  use  with  the  S  AR-based  classification  algorithm,  as  shown 
in  Fig.  3. 

It  should  be  noted  that  the  “data  points”  utilized  in  the 
curve  fitting  algorithms  were  not  simple  E-T  pairs  as  implied 
above.  Rather,  to  reduce  the  influence  of  noise  in  the  raw  E-T 
data  obtained  from  the  sensor,  and  to  limit  the  number  of  data 
points  to  a  computationally  tractable  range,  each  data  point 
was  obtained  as  an  rpproxiraation  to  the  mean  of  a  distribution 
of  individual  E-T  measurements. 

B.  Curve-Fitting  Software 

The  several  steps  outlined  in  the  above  overview  have  been 
codified  into  an  integrated  software  package  known  as  “FIT.” 
The  FIT  package  consists  of  several  dozen  individual  routines 
written  in  VAX  FORTRAN.  Included  are  routines  for 
performing  least-squares  approximations,  plotting,  error  calcu¬ 
lation,  and  generating  the  look-up  tables  required  by  the  DPU*  s 
SAR  (this  last  function  is  used  to  verify  the  performance  of 
the  on-board  table  writing  software).  The  individual  routines 
may  be  invoked  either  by  way  of  a  VAX/VMS  DCL  command 
file,  or  via  a  menu-driven  user-interface.  A  number  erf  options 
exist  for  customizing  the  general-purpose  routines  in  the  FIT 
package  to  particular  applications,  however  the  primary  utility 
of  these  routines  is  in  the  generation  of  the  coefficients  of  the 
boundary  polynomials  upon  which  the  mass  and  M/Q  group 
look-up  tables  are  based. 
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c. 


Now  define  the  column  vectors: 


The  least-squares  approximation  algorithms  utilized  in  the  FIT 
program  will  be  briefly  reviewed  here  (see  [3]  for  a  more 
extensive  development). 

Suppose  that  (xj.yj),  fa.yrf.  --  .(xN,yN)  m  independent 
measurements  (for  example,  x  =  energy  and  y  =  time-of-fligbt), 
and  that  the  distribution  is  approximated  by  the  polynomial 
y  *  P(x)  =  ci  +  cjx  + +  CjjX""!.  One  measure  of  the  total 
error  introduced  in  the  approximation  P  is  given  by  the  sum, 

N 

E(ci.c2 . c„)  =  x2  =  X  (yk  '  p<*k))2-  (1) 

k  *  1 

In  the  method  of  least-squares,  a  polynomial  P  which 
minimizes  the  error  E  is  sought  From  elementary  calculus, 
the  coefficients  of  such  a  polynomial  must  satisfy  dE/3cj  =  0 
for  i  =  1..jl  Differentiation  of  both  sides  of  Eq.  1  yields, 

N 

3E/9cj  =  -2  ]T  (yk  -  P(xk)) .  3P(xk)/3Ci 
k  =  1 


N 


=  -2 

X  (yk  -  p(*k))  • 

xkM- 

(2) 

1 

c  =  1 

N 

Hence,  3E/3cj  = 

0  if  and  only  if  V 

,(yk- 

p(xk))  xkH 

=  0. 

k  = 

1 

or  equivalently. 

N 

N 

X  ykxki'1 

=  XP(xk)* 

xki4. 

(3) 

k  =  1 

k  =  1 

which  can  also  be  written. 

N 

X  ykxkM 

= 

k  =  1 

N 

N 

N 

C1  X  xk* 

1  +c2  X  Xk‘  +  " 

•  +  Cn 

Xxkn+i-2- 

(4) 

k  =  1 

k  =  1 

1 

t  =  1 

The  final  equation  above  (Eq.  4),  taken  for  each  i  =  l...n, 
defines  a  system  of  n,  n*  order,  non-homogeneous  linear 
equations  with  respect  to  ci-c^ 

Xy*  =  CiX  1  +  ^X3*  *  -  +  CnX^”'1 
Xykxk  =  Cj£xk  +  C2Xxk2  +  "•  +  cnXxk" 


Xykxknl  “  clXxkn'1  +  'll*?  +  -  +  CnX35^'2' 


N  N  N 

Y*(  Xyk’  2>xk.  "  •  Xykxkn  l>T 

k  -  1  k  =  1  k  =  1 

and  C=*(cj.c2,„,cn)T. 

and  let  X  be  the  (n  x  n)  matrix  whose  j1^  row  is  the  vector, 

N  H  N 

(  X  xkj'1*  X  . X  *kn+j'2  )• 

k  «  1  k  =  1  k  «  1 

The  system  of  equations  given  above  can  then  be  written 
compactly  as  Y  *  XC,  and  (provided  that  X  is  invertible)  the 
solution  vector  is  simply  C  =  X_1Y.  (In  practice,  it  is  more 
economical  to  calculate  C  via  the  equation  C  =  (ZTZ)_1  Y(ZTV) 
where  Z  is  the  (N  x  n)  matrix  whose  j*  column  is  the  vector 

(xji-1,  x^-1 . “<* v  =  (yi> yi>  - .  yn)-  M  » 

easily  verified,  zTz  =  X  •\nd  ZYV  =  Y.) 

The  situation  for  multivariate  polynomials  is  entirely 
analogous  since  it  is  with  respect  to  the  coefficients  of  the 
polynomial  that  the  error  must  be  minimized.  For  example,  if 
(xj.yj.Zi).  (*2»y2»*2)*  -  .(xN.yN^N)  “>setof  measure¬ 
ments  approximated  by  a  multivariate  polynomial  of  the  form: 
z  =  P(x,y)  =  (cjj  +  C2jX  +  _.  +  Cnjx"-1)  +  (c12  +  CjjX  + 

+  CjtfX"-^  +  -.  +  (clm  +  C2rax  +  ...  +  cnmxn'1)ym*1,  then 
the  error  in  the  approximation  is  given  by, 

N 

E  =  x2  =  X  (zk  -  p(xk-yk))2-  (5) 

k  =  1 

Hence, 

N 

3E/3qj  =  -2  X  (zk  *  p(xk.yk))  •  dp(*k-ykV3cij 

k  =  1 


=  -2  X  (zk  -  p(Xk.yk))  •  Xk^ytJ’1,  C6) 

k  =  1 

N 

and  dE/Bcij  =  0  if  and  only  if  X  (zk  ‘  p(*kYk))  xki-1yic*'1 
k  =  1 

=  0,  or  equivalently, 

N  t  N 

Xvk”^1  =  X  p(xk)  *  Xk^ykJ'1-  CD 

k  =  1  .  k  =  X 


As  in  the  single  variable  case,  this  equation  leads  to  a  system 
of  n*m  non-homogeneous  linear  equations,  the  solution  of 
which  provides  the  coefficients  of  the  multivariate  polynomial 
for  which  the  approximation  error  E  is  minimal. 
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V.  Table  Writing  Software 

The  look-up  cables  required  by  DPUA's  MICS  Event  Proces¬ 
sor  are  written  by  the  SYST  able  calculation  routine.  SYST 
writes  three  distinct  tables:  The  Mass  Group  Table  used  by 
DPUA’s  Event  Processor  to  classify  events  into  mass  (M) 
groups:  the  Mass  per  Charge-State  (NVQ)  Group  Table  used  by 
the  Event  Processor  to  classify  events  into  KVQ  groups;  and 
the  Rate  Scaler  Select  Table  used  to  define  the  contents  of  the 
R0-R7  rate  scalers. 

The  procedure  followed  by  SYST  in  writing  the  Mass 
Group  Table  is  described  below  (the  M/Q  Group  Table  is 
written  in  an  analogous  fashion).  The  role  played  by  SYST  in 
constructing  the  Rate  Scaler  Select  Table  is  minimal,  and  will 
not  be  described  here. 

A.  Mass  Group  Table  Writing 

As  described  in  Section  IIB,  the  entries  in  the  Mass  Group 
Table  represent  overlayed  polynomial  curves  in  E  vs.  T  para¬ 
meter  space,  arranged  such  that  adjacent  curves  define  the  upper 
and  lower  bounds  of  a  single  mass  group.  Table  values  are 
calculated  using  coefficients  stored  in  DPU  memory.  The 
default  coefficients  are  derived  from  polynomial  fits  to  sensor 
calibration  data  as  described  in  Section  IV;  new  sets  of 
coefficients  can  also  be  uplinked  to  DPUA  using  ground 
commands. 

The  Mass  Group  Table  is  used  by  the  hardware-based 
Successive  Approximation  Routine  of  the  Event  Processor  to 
classify  each  event  into  one  of  up  to  16  mass  groups,  based  on 
the  8-bit  energy  (E)  and  10-bit  time-of-flighi  (T)  values 
obtained  from  the  sensor.  The  table  is  divided  into  16  sections, 
one  per  mass  group,  each  of  which  contains  256 10-bit  entries. 

Each  mass  group  section  of  the  table  is  indexed  by  energy 
values  ranging  from  E  =  OtoE  =  255.  For  each  energy  value 
E,  a  corresponding  nominal  time-of-flight  value  T  is  calculated 
by  SYST  as  T  =  2l,  where  t  =  P(x)  =  co  +  cjx  +  C2X2  +  tgx3  Fig.  4.  Mass  Group  Classification  Schemes, 

and  x  =  Log2E,  and  then  the  complement  of  this  value  relative 
to  1024  is  written  to  the  Mass  Group  Table.  A  (Efferent  set  of 

coefficients  is  used  to  calculate  the  table  values  for  each  mass  fication  scheme  in  which  all  events  with  sufficiently  small 
group.  The  coefficients  define  non-intersecting  curves  in  E  vs.  energy  values  are  placed  in  “special”  mass  groups.  Either  two 
T  parameter  space  (see  Fig.  3)  which  represent  the  boundaries  or  three  special  mass  groups  may  be  allocated.  With  two 
between  successive  mass  groups  (fifteen  such  curves  are  special  mass  groups  (Fig.  4b),  every  event  with  energy  less 
therefore  necessary  to  partition  the'space  into  sixteen  mass  than  or  equal  to  the  pre-defined  constant  ECUT  is  placed  in 
groups).  In  particular,  the  time-of-flight  values  in  a  given  either  Mass  Group  0  or  Mass  Group  I,  depending  on  whether 
mass  group  section  of  the  table  represent  the  lower  boundary  the  corresponding  time-of-flight  value  is  less  than  or  equal  to 
of  the  mass  group,  hence  no  values  need  be  written  for  Mass  the  constant  TCUT1  (Mass  Group  0),  or  greater  than  TCUT1 
Group  0.  (Mass  Group  l).  Events  with  energy  values  greater  than 

In  the  basic  classification  scheme  up  to  16  mass  groups  ECUT  are  classified  as  in  the  baric  scheme  described  above.  In 
may  be  defined  via  stored  coefficients  as  above.  It  is  important  particular,  if  an  event  with  E  >  ECUT  falls  into  either  Mass 
to  notice  that  fewer  than  16  mass  groups  may  be  allocated,  Group  0  or  Mass  Group  1  under  the  baric  scheme,  then  it  will 
simply  by  specifying  appropriate  coefficients.  For  example,  also  be  placed  there  under  the  special  scheme,  thereby  cor- 
specifying  “barrier-down  coefficients”  for  Mass  Group  12  rupting  the  interpretation  of  these  mass  groups  as  low-energy 
ensures  that  only  the  lowest  12  mass  groups  will  be  used,  event  collectors.  To  avoid  this  problem,  the  coefficients  for 
These  coefficients  result  in  zeros  being  written  to  the  Mass  the  Mass  Group  2  section  of  the  table  should  be  changed  to 
Group  Table  for  the  specified  group;  hence,  since  for  any  T,  “barrier-up  coefficients”.  These  coefficients  cause  a  saturated  T 
T  +  0  5  1023,  no  event  will  ever  be  placed  in  Mass  Group  12  value  to  be  written  to  the  table  locations  corresponding  to 
or  any  higher  numbered  group,  i.e.,  only  Mass  Groups  0  E  >  ECUT  in  the  Mass  Group  2  section  of  the  table,  the  net 
through  1 1  will  be  used.  result  of  which  is  that  all  events  with  E  >  ECUT  which  would 

In  addition  to  the  basic  classification  scheme  outlined  above  have  fallen  into  Mass  Group  0  or  Mass  Group  1  under  the 
and  illustrated  in  Fig.  4a,  SYST  also  offers  a  special  classi-  basic  classification  scheme  are  instead  placed  in  Mass  Group  2. 
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The  situation  with  three  special  mass  groups  is  similar 
(Fig.  4c)  except  that  an  additional  constant,  TCUT2,  is  used  in 
conjunction  with  TCUT1  to  place  events  with  E  <  ECUT  in 
either  Mass  Group  0  (TS  TCUT1),  Mass  Group  1  (TCUT1  < 
T  <  TCUT2),  or  Mass  Group  2  (T  >  TCUT2).  As  noted 
above,  to  maintain  the  integrity  of  the  3  special  mass  groups, 
the  coefficients  for  both  Mass  Group  2  and  Mass  Group  3 
should  be  changed  to  “barrier-up  coefficients”. 

SYST  also  offers  an  option  through  which  all  events  with 
energy  equal  to  the  maximum  value  of  25S  are  placed  in  Mass 
Group  13  (shown  in  Fig.  4d  together  with  the  special  classi¬ 
fication  scheme).  This  option  may  be  used  in  conjunction 
with  any  of  the  above  classification  schemes,  that  is, 
independent  of  the  number  of  special  mass  groups,  if  any. 
However,  it  should  be  noted  that  if  the  selected  classification 
scheme  makes  use  of  all  16  mass  groups  (including  Mass 
Group  15),  then  Mass  Group  15  will  collect  lower  energy 
events  in  addition  to  events  with  E  =  255.  Thus  this  option  is 
only  used  when  “barrier-down  coefficients”  are  simultaneously 
being  employed  to  deflect  lower-energy  events  away  from 
Mass  Group  15. 

VI.  Summary 

The  VIKING  MICS  DPU,  the  CRRES  LOMICS  and  MICS 
DPUs,  and  the  forthcoming  CAMMICE  DPUs  all  employ  the 
M  vs.  M/Q  analysis  techniques  described  here.  The  combined 
use  of  M  vs.  M/Q  scalers,  rate  scalers,  and  sample  direct 
events  allows  event  parameter  space  to  be  monitored  at  a  range 
of  parametric  and  temporal  resolutions  within  the  constraints 
imposed  by  the  limited-bandwidth  telemetry  channel  of  the 
spacecraft.  Scaler  compression  further  conserves  telemetry 
bandwidth,  without  significant  loss  of  precision. 

A  hardware-implemented  binary  search  strategy  called  the 
Successive  Approximation  Routine  (SAR)  is  used  to  classify 
events  into  mass  and  M/Q  groups  using  memory-resident 
tables  written  by  the  DPU  software.  Use  of  the  SAR-based 
scheme  results  in  a  significant  reduction  in  memory  require¬ 
ments  compared  with  direct  table  look-up  techniques,  without 
creating  a  processing  bottleneck  at  the  DPU. 

The  tables  used  by  the  SAR  are  user-configurable.  A  num¬ 
ber  of  options  exist  whereby  events  with  significantly  small  or 
large  parameter  values  may  be  placed  in  “special”  groups,  seg¬ 
regated  from  the  general  event  population.  The  tables  may  be 
reconfigured  in-flight  by  simply  changing  the  variables  con¬ 
trolling  operation  of  the  SYST  table  writing  routine. 

The  discrete  polynomial  approximations  contained  within 
the  tables  were  obtained  using  a  novel  combination  of  local 
(single  group)  and  global  curve-fitting  algorithms.  The 
resulting  approximations  are  globally  consistent  in  the  fashion 
required  by  the  SAR-based  classification  scheme. 


VII.  Conclusion 

Clearly  techniques  similar  to  those  described  here  are  necessary 
to  reduce  the  enormous  quantity  of  sensor-generated  data  and 
maximize  the  utilization  of  the  available  telemetry  bandwidth. 
Although  other  on-board  analysis  and  reduction  schemes  might 
be  considered,  the  M  vs.  M/Q  classification  technique  utilized 
by  the  VIKING  and  CRRES  DPUs  has  proven  to  be  quite 
successful,  and  reasonably  simple  to  implement.  Slight  modi¬ 
fications  to  the  existing  algorithms  are  anticipated  for  the 
forthcoming  CAMMICE  DPUs,  including  the  addition  of 
several  new  options  to  the  SYST  table  writing  routine  to 
allow  further  specialization  of  mass  and  M/Q  groups.  The 
curve-fitting  routines  are  also  being  refined  to  provide 
additional  tools  for  evaluating  the  “goodness"  of  the  resultant 
classification  scheme.  In  general,  the  proposed  modifications 
seek  to  add  greater  flexibility  to  the  existing  algorithms, 
without  altering  the  essence  of  these  flight-tested  techniques. 
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TECHNOLOGY  OPERATIONS 


Tbe  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  Tbe  Corporation's  Technology 
Operations  supports  tbe  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of 
new  materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new 
forms  of  carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment;  contamination,  thermal  and  structural  control;  high 
temperature  thennomechanics,  gas  kinetics  and  radiation;  lubrication  and  surface 
phenomena. 

Space  and  Environment  Technology  Center:  Magnetospberic,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


